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Abstract 

Usual introductions of the concept of motion are not well adapted 
to a subsequent, strictly tensorial, theory of elasticity. The consider- 
ation of arbitrary coordinate systems for the representation of both, 
the points in the laboratory, and the material points (comoving co- 
ordinates), allows to develop a simple, old fashioned theory, where 
only measurable quantities — like the Cauchy stress — need be intro- 
duced. The theory accounts for the possibility of asymmetric stress 
(Cosserat elastic media), but, contrary to usual developments of the 
theory, the basic variable is not a micro-rotation, but the more fun- 
damental micro-rotation velocity. The deformation tensor here intro- 
duced is the proper tensorial equivalent of the poorly defined defor- 
mation "tensors" of the usual theory. It is related to the deformation 
velocity tensor via the matricant. The strain is the logarithm of the de- 
formation tensor. As the theory accounts for general Cosserat media, 
the strain is not necessarily symmetric. Hooke's law can be properly 
introduced in the material coordinates (as the stiffness is a function 
of the material point). A particularity of the theory is that the compo- 
nents of the stiffness tensor in the material (comoving) coordinates are 
not time-dependent. The configuration space is identified to the part 
of the Lie group GL + (3) that is geodesically connected to the origin 
of the group. 
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1 Introduction 



Elasticity has often been the model theory for building other theories. For 
instance, Maxwell used the elastic analogy to develop his equations, and 
Elie Cartan was inspired by elasticity 1 to propose the Einstein-Cartan grav- 
itation theory 2 . 

Today, elasticity is mainly being developed by applied mathematicians, 
whose principal goal is mathematical generality and consistency, at the cost 
— in my opinion — of blurring the difference between the tensors beloved 
by Cauchy and Einstein 3 , and their generalization as abstract mappings 
between linear spaces (subjected to pull backs and pull forwards). Also, 
some confusion exists on the physical interpretation of the the different 
stress "tensors" introduced in the usual theory, that are better seen as just 
auxiliary computational tools. Therefore, in this note I take the notion of 
tensor much more basically than usual texts in deformation theory. And 
I resolutely take the old-fashioned index-based notation: it is always easy, 
once the basic mathematics and physics are understood, to move towards 
more abstract notations. 

Because I want to highlight some simple messages, I completely ignore 
the dynamic part of the problem (conservation of mass, of linear momen- 
tum, and of angular momentum are not considered): only the definition of 
strain, of stress, and of internal elastic energy are examined. And, of course, 
the relation between stress and strain. 

We are now celebrating the one-hundredth anniversary of the ground- 
breaking work of the Cosserat brothers (Cosserat and Cosserat, 1909). Most 
of the scientists familiar with the Cosserat theory understand that the usual 
assumption of symmetric stress breaks the beauty of the mathematics, in 
exactly the same way as anyone who understands the Einstein-Cartan the- 

1 Cosserat theory, with asymmetric stress. 

2 While in Einstein's theory space-time may have curvature but no torsion, in the Eistein- 
Cartan theory, both curvature and torsion may exist, torsion allowing to take into account 
the existence of spin in realistic models of matter. 

3 A tensor is an intrinsically defined quantity at some point of the differentiable manifold 
representing the space-time. It belongs to one of the linear spaces that can be built by dif- 
ferent tensor products of the tangent linear manifold and its dual. The sum and the product 
by a real number are the two basic operations for tensors. In reality, there are objects — like 
a rotation operator — that share with tensors the property of being intrinsically defined, but 
whose natural operations are the product and the exponentiation to a real number. These 
are not — strictly speaking — tensors, and are better seen as points on some Lie group tan- 
gent to the space-time manifold (Tarantola, 2006). Their logarithm is almost a tensor, as the 
sum and product for a real number is defined, but this sum is not commutative. So the 
reader may now understand how restrictive this author is in the use of the term tensor. 
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ory of gravitation knows that assuming vanishing space-time torsion (thus, 
a symmetric connection) takes out most of the beauty of the Bianchi con- 
servation equations. 

So, in this theory, the stress tensor can be asymmetric. And, following 
the Cosserats, we interpret the antisymmetric part of the stress as the effect 
of micro-rotations (of, say, the material "molecules"). Yet, I choose not to 
follow common practice of explicitly considering the micro-rotations as a 
primitive ingredient of the theory. For a rotation is always relative to some 
initial configuration, whereas a rotation velocity is not. And, for the same 
reason, (symmetric) strain is not a primitive ingredient; the deformation ve- 
locity is. This matters, as the deformation velocity is not simply the strain 
rate. 

2 Movement 

2.1 Frame of reference 

Let us start by assuming the existence of Galilean frames of reference, and 
by choosing a particular one, say Q , with respect to which all tensor fields 
(velocities, stresses, etc.) are defined. The time coordinate t is Newtonian 
time. We do not need to assume that the spatial part of the Galilean frame Q 
is Euclidean 4 , or that it is necessarily thee-dimensional (although I shall use 
a language adapted to the three-dimensional case). A space point of Q may 
be denoted using a letter like T . 

A tensor field may be represented using a notation like 

s = <r(0>,0 . (1) 

Here, s denotes the tensor at space- time point {7, t} while a denotes 
the function of the space-time coordinates, exactly as when — in elementary 
mathematics — one writes y = f(x) . 

2.2 Motion 

A deforming body B is assumed to occupy the whole 5 of the space. Its 
points, called material points, are assumed to be individually identifiable 

4 What is convenient if the theory is to be applied to some non-flat submanifold of the 
three-dimensional physical space. 

5 The modifications to be made when the body B occupies only part of the space are 
quite trivial conceptually although technically complex. 
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and their trajectory in the Galilean frame Q knowable (at least, in principle). 
The trajectories of the material points is the function 

y = /(M,t) , (2) 

specifying, at every instant t , the laboratory position 7 of any material 
point M . The continuity hypothesis is that the inverse function exist: 

M = F(3>,f) . (3) 

With this, equation (1) can be completed by introducing a new function: 

s = a(y,t) (9 = f(M,t)) 

= E(M,f) (M = F(:P,f)) • 



2.3 Laboratory coordinates 

The theory could be developed using intrinsic notations only, i.e., without 
writing tensor equations in terms of the components of the tensors in the 
natural basis associated to some coordinates. But it is well-known that, 
as far as the coordinates are arbitrary, component-based tensor equations 
are intrinsic. Many mathematicians prefer more abstract, component-free, 
notations, and there is no problem with that, excepted that pedagogy may 
command using expressions that are as explicit as possible. Sometimes, one 
fails to understand that coordinates are a tool for discovery: there are many 
examples where it is the discovery of a system of coordinates adapted to a 
problem that has allowed its resolution 6 . 

So, let us assume that some (fixed, arbitrary) coordinate system, x = 
{x 1 } = \x x ,x 2 , x 3 } is chosen in the spatial part of the Galilean frame of ref- 
erence 7 . This coordinate system, together with the Newtonian time t con- 
stitute what we shall call a system of (space-time) laboratory coordinates. At 
every space point x the natural vector basis e,(x) is considered, together 
with the associated natural tensor basis e;(x) <S> e,-(x) ...e k (x)<S> e e (x) . . . . 
Then, for the tensor field in equation (1) we can now write: 

s = s{x,t) = SijJ e -(x,t) e'(x) ®e>(x) ...e k (x) ®e e {x) ... . (5) 

6 The discovery of the Schwarzschild coordinates in space-time was fundamental for the 
discovery of spherically symmetrical solutions to Einstein's equations. More recently, it was 
the discovery of a space-time coordinate system that allowed to properly formulate — and 
solve — the problem of a fully relativistic positioning system (Coll, 2002; Tarantola et al., 
2009). Not to speak about Lagrange's discovery of the material coordinates. . . 

7 The space is not assumed to be Euclidean, so, a fortiori, these coordinates are not as- 
sumed to be Cartesian. 



5 



Note that, by definition of the laboratory coordinates, the vector basis e, (x) 
is time-independent (contrary to the comoving vector basis about to be in- 
troduced). 



2.4 Material (comoving) coordinates 

Let us now assume that an arbitrary system of material coordinates X — 
{X 1 } = {X 1 , X 2 , X 3 } has been chosen, as suggested in figure 1. By defini- 
tion, the material coordinates of any material point M have constant val- 
ues. The trajectories in equation (2) have now the more concrete, coordinate- 
based, representation x i = ^'(X^X 2 , X 3 , t) , (i = 1,2,3) , or, for short, 

x = <p(X,t) . (6) 

The continuity hypothesis is that these functions are continuous and invert- 
ible, so that the inverse functions exist: X 1 = x 2 ,x 3 , t) , (J = 1,2,3). 

For short, we simply write 

X = <D(x,f) . (7) 

Seen from the laboratory (i.e., from the Galilean frame of reference), the 
material coordinates deform. Therefore, the natural vector basis associated 
to the material coordinates (these vectors — as all tensors of the theory — 
are defined with respect to Q ) is time-dependent, so the notation e/(X, t) 
has to be used. The availability of the material system of coordinates, and 
of the associated tensor basis, allows to further complete equation (5), in- 
troducing also the material components of the tensor field: 



s(x,t) 

Siij'-(x,t ) e\x) ® e'(x) . . . e k (x) ® e e (x) . . . 



= S(X,f) 



(8) 



= s 



I]... 



KL... 



X, t) e 1 (X, t) ® e 7 (X, t ) . . . ej<(X, t) ® e L (X, f ) 



with the understanding that for these identities to make sense one has to 
use the replacements in equations (6)-(7). 

2.5 Coordinate change 

Associated to the coordinate changes in equations (6)-(7) are the matrices 
of coefficients 

A ! i(X,0 = Hj(X,0 ; \ l i{x,t) = ^r(^0 / ( 9 ) 
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Figure 1: The "laboratory" is assumed to be a Galilean frame of reference. 
Two snapshots of the laboratory, at two instants t\ and f 2 / with a deform- 
ing body occupying part of the laboratory space. In red, a system of labora- 
tory coordinates, x = {x 1 } , and, in blue, a material (or comoving) system 
of coordinates, X = {x 1 } . The physical space is not necessarily Euclidean, 
and both coordinate systems are arbitrary. The motion of the body is rep- 
resented by the functions x = cp(X,t) , prescribing, for any time t, the 
laboratory position x of every material point point X . This author be- 
lieves that the usual drawings, with plenty of axes and vectors, involve 
superfluous notions. 
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that are mutually inverse: A'k A K j = <5'j , A\ A j = 8 1 ] . A well-known 
result from tensor calculus is that the vector and form bases in equation (8) 
are related as 8 ej = A ! je,, e 1 = A 1 , e' r e ; = A^ e/, e ! = A'/e I / from 
where follows that the components Sij„ ke -(x,t) are related to the compo- 
nents s J/ .. KL -(X,f) as 9 

Sl j. KL - = X^tij ...A K k \ L i ...s ij J i ~ , (10) 

or, equivalently as 10 

Sy... = A A'j . . . A K A L ■ ■ ■ SlJ... ■ (11) 



2.6 Metric 

We shall assume that the components of the metric gij{x) are known in 
the laboratory coordinates. Quite often, the space is going to be Euclidean, 
and, in this case, the gij(x) are just the components of the Euclidean metric 
in the (arbitrary) laboratory coordinates, but let us not assume that we are 
in this special situation. We can naturally write 

g = gW = gij(x) e i {x)^e'{x) , (12) 

denoting by the symbol g the metric tensor itself. Note that the compo- 
nents gij(x) of the metric tensor in the laboratory coordinates are time- 
independent. This is not so in the material coordinates, where one has 

g = G(X) = gI j(X,t)e\X,t)®e?(X,t) , (13) 

as both, the basis ej(X, t) and the components gij{X, t) , are time-depen- 
dent. The components gjj can be expressed as a function of the compo- 
nents gij as gij = A 1 i A 1 j gjj, i.e., more explicitly, 

gIJ (X,t) = A i I (X,t)Aij(X,t)g lj ( ( j ) (X,t)) . (14) 

8 I.e., more explicitly ej(X,f) = A* x (X,f) e ; (^(X,f) ) , e z (X,i) = A 1 ,{cp(X,t) , t ) 
e'(0(X,i)), e ; (x) =A I i {x)e I (<S>{x,t),t), e'(x) = A'j( ®(x, t) , t) e 7 ( 3>(*, t) , t). 

9 I.e v more explicitly, s;;.. KL ~(X,t) = A'j(X, t) A?/(X, t) . . . A K k ( <5(x, t) , t) A L e { 0>(x, 

0,0 ... Slj .. ke -(<i>(x,t),t). 

10 I.e., more explicitly SijJ d -(x,t) = A 1 ,{x,t) A^{x,t) . . . A k K ( ®(x, t) , t) A f L ( ®(x, f) , 
... Sl j.. KL -(^(x,t),t). 
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2.7 Velocity 



The considered motion defines a velocity field, namely, the velocity of all 
the material points with respect to the Galilean frame of reference, 

V 1 = dft/dt , (15) 

an equation that maybe better understood if all the variables are explicited: 
v'(x, t) = (dcp'/dt) ( 0(x, t) , t) . To obtain the expression of the velocity 
field in the material coordinates, one can just use v 1 = A ,• v l . Alterna- 
tively, one has 

v l = -W/dt , (16) 

i.e., v\X, t) = — (90V dt) ( <p(X,t) , t) . The covariant components of the 
vector field are Vi(x, t) = gij(x) vi(x,t) and Vj(X,t) = gjj(X,t) v^(X,t) . 



2.8 Deformation velocity 

The notion of "strain tensor" is subtle, and only to be introduced later. 
A robust notion is that of deformation velocity (similar, but different from 
the "strain rate" to be later introduced). The deformation velocity tensor, 
denoted d , can be introduced using any of the two equivalent definitions 

dij = I {Vivj + Vjvj) ; dij = \ (V(Dj + Vpi) . (17) 

Because of the definition of material coordinates, one has the property 11 



d u (X,t) = l gIJ (x,t) 

The vorticity 



(18) 



cd u = \ (y lVj - V/Pj) ; CDij = \ {ViVj-VjVi) , (19) 

represents a local "mesoscopic rotation velocity". It has not to be mis- 
taken for the fundamental "microscopic rotation velocity", about to be in- 
troduced. 



11 This can be obtained by evaluating the partial time derivative of expression (13). As 
3fG = 0, the result follows when using the property 3fe 7 = — (Vx» I )e K . 



9 



2.9 Rotation velocity 

So far, the considered motion has only considered the "translational" move- 
ments of the material points, considered as featureless. More realistic con- 
tinuous models of matter also consider the possibility that the individual 
material points (i.e., the "molecules") can rotate. This is, for example, the 
case for the fluids where the existence of a spin density is to be consid- 
ered 12 . It is also the case in the theory of elastic media where the stress 
tensor is not assumed to be symmetric. This theory, developed by Eugene 
and Francois Cosserat (Cosserat and Cosserat, 1909), also contains micro- 
rotations. 

My goal in this note is the study of general elastic media, not of fluids 
with spin. But I think it is a mistake to develop the theory of asymmetric 
elasticity starting with the notion of rotation. For the notion of (instanta- 
neous) rotation velocity is more primitive. If necessary, the rotations have to 
be evaluated by properly integrating the rotation velocity. 

So, let us consider (Cosserat media) that the material points (or "mo- 
lecules"), in addition to their translational motion, may rotate, i.e., every 
material point shall have associated a rotation velocity. This corresponds 
to an antisymmetric tensor field: 

co = Q(X,f) = co IJ (X,t) e'iXj) ^)e J (X,t) 
= w{x,t) = C0ij(x,t) e ! (x) <g> e'(x) , 

with co j j + cojj = and cOjj + cou = . Again, this intrinsic (or "micro- 
scopic) rotation is not to be mistaken for the vorticity (equation (19)). 

2.10 Movement velocity 

The deformation velocity tensor d is, by definition, symmetric. The sum 
of the symmetric deformation velocity and of the antisymmetric rotation 
velocity, 

A;/ = dij + coij ; Aij = dij + cojj , (21) 

is going to be called the movement velocity. This tensor A has no special 
symmetry. 

12 For the theoretical beauty of a relativistic theory of fluids with spin, see Halwachs 
(1960). 
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3 Stress 



In this text, the symbol s represents the Cauchy stress tensor, first intro- 
duced by Augustin Cauchy around 1822 (Cauchy 1841) sometimes called 
the physical stress tensor. Its components are defined via 

s = S(X,f) = s I J (X,0e I (X,0®e^(X / 
= s(x,t) = s'j(x,t)ei(x)(3ei(x) . 

In a theory like this one, where the term tensor is used in a very restrictive 
sense, other matrices of numbers have no place, as, for instance, the differ- 
ent Piola-Kirchhoff stress "tensors" (Truesdell and Toupin, 1960; Eringen, 
1962; Malvern, 1969; Marsden and Hughes, 1983). 

In the general theory developed here is is not assumed that the stress 
tensor is symmetric. So, generally, Sjj ^ Sji , and s» ^ su . 



4 Viscosity 

Having introduced the stress tensor s and the movement velocity ten- 
sor A , we can introduce the notion of linear viscosity by just assuming 
proportionality between the two: 

s = E : A . (23) 

Here, E is the viscosity tensor, a positive definite tensor with some symme- 
tries 13 . I am not going to further develop here the linear theory of viscosity. 



5 Elasticity 

5.1 Dealing with two temporal variables 

Some of the "tensor functions" to be introduced below are defined with re- 
spect to some reference time, that we will denote to . To denote such tensor 
functions, we shall use the notation s(X, t; to) (note the " ; "), this meaning 
that to is considered to be a fixed constant. In particular, no time deriva- 
tives can be considered with respect to to . As in material coordinates, the 
natural basis is time dependent, it will always be considered that the vector 
basis to be used is that at t . For instance, the component of a tensor t as 

13 The first group of two indices and the second group of two indices can be permuted. 
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the deformation or the strain tensor are t'j(x, f; to) in the laboratory coordi- 
nates and t I j(X, t; to) in the material coordinates, in the precise sense that 
one has 

t = f;*o) e;(x) ®e>(x) x = <p(X,t) 

= t'f(X,f;to)e f (X,t)®e'(X,t) X = &(x,t) . 

With this convention, objects like the deformation tensor or the strain ten- 
sor (that depend on the parameter to ) are ordinary tensors. In particular, 
from the relations (space variables omitted) 

e, = A^OeiCO ; e,(f) = A ! ^)e ; - , (25) 

the usual rule for the change of components under a change of coordinates 
follows (space variables omitted): 

tJ(t;to) = K l i{t)t^{t;to) \1j(t) 

j T (261 

ti'{t;to) = A',(() t/(f;i ) A J j(t) . 

5.2 Deformation 

I shall now introduce a tensor denoted using the symbol F . Its components 
in the laboratory and the material coordinates are defined via 

F = Fj(x,t)e i (x)®e j (x) x = <p{X,t) 

= Fi ] {X,t) e\X,t) <8>e 7 (X,f) X = 0(x,f) . 

There are two equivalent definitions, one using the laboratory coordinates, 
and one using the material coordinates (space variables omitted): 

Fj(t;t Q ) = \> K (t) \ K i(t ) ; F/(f;f ) = ^i(t) A^(t ) . (28) 

These two definitions are equivalent in the sense that they satisfy the tensor 
rule expressed by equation (26). 

The tensor F so introduced is the proper tensor replacement for the de- 
formation gradient "tensor" of the usual theory. We shall acall F the defor- 
mation gradient tensor (without the quotation marks). There is no risk of 
confusion with the common deformation gradient, as, while that object has 
mixed indices, like in F; ' , our deformation gradient tensor always has in- 
dices Fjl (in laboratory coordinates) or Fjl (in material coordinates). 

We now need to make a digression. While sometimes the symbol "trans- 
pose" is used by analogy with matrix theory, we need here to be precise. In 
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particular, we need to carefully define the adjoint *t of a tensor t, as this 
is done in appendix 11.2. Applying this general definition to our present 
problem, where we have two coordinate systems — the laboratory one and 
the material one — we find the two relations 

*f f ;(f;fo) = g ik tk e (t;t )gej ; *f 7 /(f;fo) = g IK (t) t K L (t;t ) gLJ (t) . (29) 

These two definitions are equivalent in the sense that they satisfy the tensor 
rule expressed by equation (26). 

We can now introduce a fundamental tensor of deformation theory, that 
we shall call the squared deformation tensor. 

C = *F F . (30) 

Explicitly, this is 

Cj(t;to) = gikF/(t;t Q )g fm F m j(t;t ) 

Cj(t'.k) = giK(t) F L K (t;t ) g LM (t) F M J(t;t ) ^ 

These two definitions are equivalent in the sense that they satisfy the ten- 
sor rule expressed by equation (26). This tensor C satisfies the following 
properties: 

Property #1: The squared deformation tensor is self-adjoint, i.e., one has 

*C = C . (32) 

(I leave to the reader to express this equation in both, laboratory and mate- 
rial coordinates.) 

Property #2: The squared deformation tensor is symmetric, i.e., when defin- 
ing 

Cij{t;t ) = Ci k {t;t ) g kj ; C I} (t;t ) = C/(f;f ) g KJ (t) , (33) 
one has 

Cij(t;t ) = Cji(t;t ) ; Cij(t;t ) = C n {t;t ) . (34) 



Property #3: In material coordinates, the components of the squared defor- 
mation tensor can be expressed as (making explicit the space variable X ) 



C/(X,f;f ) = gi K (X,t)g K J(X,t Q ) 



(35) 
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In the laboratory coordinates, no special simplification occurs, so one just has 

cHt;t ) = A',-(0 Ci'(f; fo)A>/(0 . (36) 

These three properties are easily demonstrated, via direct substitution. 

While the components of the squared deformation tensor C in the ma- 
terial coordinates, C/ , correspond to the usual definition of the right Cau- 
chy-Green deformation "tensor", the components in the laboratory coordi- 
nates, Cji , correspond to the usual definition of the left Cauchy-Green de- 
formation "tensor". While in the conventional theory two different names 
(right- and left- deformation tensor) are used, as well as two different sym- 
bols (usually C and B ), we see that, in reality, there is only one tensor 
(with, of course, different components in different bases). The deformation 
tensor originally introduced by Cauchy (in 1828) is, in fact, the inverse of 
our squared deformation tensor C . 

In a work like this one, it is out of question to give different names to 
a unique tensor, so we have to use a single name for the tensor C . As in a 
one-dimensional elongation problem, the determinant of this tensor is 

detC = (£(t)/£(t Q )) 2 , (37) 

it seems that the name here used (squared deformation tensor) is adequate. 

Property #4: Using (18) and (35), one immediately obtains 

C I K {X,t;t )C K J (X,t;t ) = 2d/(X,f) , (38) 

where a "dot" denotes partial time derivative, and where the denote the com- 
ponents of the tensor C 1 . 



The square root of the squared deformation tensor, 



D = VC 



(39) 



shall naturally be named the deformation tensor. Is is also symmetric and 
self-adjoint. 

Property #5: From equation (38) it follows (using the fact that d and C are 
symmetric) 



D I K {X / f r t )D K J{X,t;t ) = dJ{X,t) . (40) 
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The expression in equation (40) mathematically corresponds to the no- 
tion of declinative (Tarantola, 2006), that is the proper time derivative to be 
introduced for this kind of tensors 14 : the declinative of the deformation is the 
deformation velocity. 

Property #6: From relation (40) it follows, using the matricant theory (see 
appendix 11.1), that one has (variable X implicit) 

t t t' 

D/(t;t ) = V + [dt'd I J{t')+ Ut'di K {t') f dt" d K ] \t") + . . . . (41) 

J f o J t o J to 

So, while property #5 says that the deformation velocity is a "properly 
defined" time derivative of the deformation tensor, this property #6 gives 
the inverse relation, expressing the deformation tensor as a "properly de- 
fined" time integral of the deformation velocity tensor. In section 8 we 
shall see that this is an integration of a Lie group manifold (representing 
the configuration space), with continuous parallel transport to the origin of 
the group. 

Properties #5 and #6, taken together, suggest that our deformation ten- 
sor D is intimately connected to the deformation velocity tensor d . The 
deformation tensor D is, therefore, a fundamental tensor in the theory of 
continuous media. 

Property #7: One has 15 

detD(X, f;f ) = exp / dt' trace d(X, t') . (42) 

Jt 

As det D expresses the ratio between final and initial volumes, this relation 
relates that ratio to the time integral to the trace of the deformation velocity 
tensor. 

5.3 Symmetric strain 

Cauchy originally defined the strain as 

E = |(C- I) , (43) 

14 I am reluctantly using the name tensor here (see footnote 3). 
15 See appendix 11.1. 
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but many lines of thought suggest that this was just a guess, that, in reality, 
is just the first order approximation to the more proper definition 

E = log^C=i(C-l)-i(C-l) 2 + ... , (44) 

i.e., in reality, 

E = log D = ( D - I ) - i( D - I ) 2 + . . . . (45) 

But this requires some care, as the logarithm of a real matrix is not always 
real. 

Definition (symmetric strain): Let be D a (symmetric) tensor 16 that be- 
longs to the part of the Lie group manifold GL + (3) that is geodesically connected 
to the origin of the group. Then (Tarantola, 2006), the logarithm of D is a real 
tensor 17 , and the (symmetric) strain associated to D is defined as 

E = log D . (46) 



The reason for the strain being not defined for an arbitrary D is ex- 
plained in section 8. The strain defined logarithmically is often named 
natural strain of Hencky strain (e.g., Truesdell and Toupin (1960), Rougee 
(1997)). 

The components of the strain tensor are, of course, defined via 

E = Ej(x,t;t ) e i (x)®e j (x) x = <p{X,t) 

= E I J (X,t;t )e I (X,t)®ej(X,t) X = <f>(x,t) . 1 ° 

This is a bona-fide tensor. It is easy to see that this strain tensor is both, 
symmetric and self-adjoint. 

An actual computation of the symmetric strain can be done in both, the 
laboratory and the material coordinates. First, one may use the property, 

E = logv^C = UogC , (48) 

so one does not have to care about the square-root. Then, computing the 
logarithm of a tensor just amounts to compute the logarithm of a matrix 

16 Or, if the reader prefers, the matrix representing the covariant-contravatiant compo- 
nents of the tensor in some basis. 

17 I.e., the matrix with the components is real. 
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whose entries are the mixed components (i.e., covariant-contravariant or 
contravariant-covariant) of the tensor in any basis. The result so obtained 
is intrinsic (i.e., independent from the basis being used) 18 . 

There are different ways for computing the logarithm of a second-rank 
tensor given its mixed components. These range from the series expan- 
sion 19 to the fully analytical methods proposed in Tarantola (2006). 



5.4 Asymmetric strain 

In equation (21) I have introduced the movement velocity A as the sum 
of the (symmetric) deformation velocity and the (antisymmetric) rotation 
velocity: 

A;/ = dij + coy ; A;y = dij + wij . (49) 

To introduce the notion of an asymmetric strain, we just need to col- 
lect some of the equation above, and drop the assumption that tensors are 
symmetric. 

The symmetric deformation tensor D generalizes into the asymetric 
deformation tensor A that bears with A , the same relation that D bears 
with d . The equivalent of equation (40) is 



A I K (X / t;t )A K I(X,t;to) = Ai J (X / t) , (50) 



while the equivalent of the relation (41) is 

rf 

c 

'to 



Aj(t;t ) = Si J + [dt'Aj(t')+ fdt'Afit') fdt"& K l(t") + ... . (51) 



The components of the tensor A in the laboratory coordinates are to be 
obtained via the usual relation implied by a change of coordinates: A$ = 

The equation defining the asymmetric strain is just the equivalent of 
equation (46): 



E = log A . (52) 

We could use a different symbol for the asymmetric strain, but as this is just 
an "obvious" generalization, let us keep the same symbol E . As above, the 

18 This follows directly from the property that, for any invertible matrix M , and for any 
matrix C, one has M (logC) IVT 1 = log(M C M" 1 ) . 

19 It may be that the expansion (logC)/ = (C„ 6 - S a b ) - \ (C a c - S a °) (C c b - 8 c b ) + 
j (C fl c — 5 a c ) (C c d — 5 c d ) (Q fc — 8f) — . . . converges. This, of course, is nothing but 
logC = ( C -I)-i(C-I) 2 + i(C-I) 3 -.... 
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strain is only defined if A belongs to the part of the Lie group manifold 
GL + (3) that is geodesically connected to the origin of the group, i.e., in 
fact, if log A is real. 

In the situation where there are no micro-rotations, co = , so A = d , 
and A is symmetric. This is obviously the special case analyzed in sec- 
tion 5.3 (symmetric strain), and we do not need to return to it. 

Let us then analyze the other extreme situation, where there are only 
micro-rotations. Then, the (symmetric) deformation velocity d is zero, 
A = co , and A is antisymmetric. The matricant series (51) then gives 
RiJ (X, f; to) the (orthogonal) operator representing the total rotation 20 be- 
tween fo and f: 

t t t' 

R/(t;t„) = V+ fdt'tofit') f dt" tu K J(t") + ... . (53) 

Jtg J to J tg 

Note how the usual Cosserat micro-rotation enters the scene in the theory 
here proposed: as a (quite complex) quantity derived — via the matricant — 
from the (more elementary) micro-rotation velocity. 



5.5 Hooke's law 

During the evolution of a deforming medium, different values of the time f 
are considered. In elasticity, one assumes that there is some reference "con- 
figuration" that is kept in memory by the deforming medium. Let us as- 
sume that this is the configuration at instant fo , and let us simplify the 
theory by assuming that there is no "pre-stress", i.e., that the stress at in- 
stant fo is zero. To remember that special condition, let us, from now on, 
change our notation for the stress, and denote it s(f; fo) instead of just s(f ) . 
The initial condition is then 



s J /(f ;fo) = 



(54) 



The proper formulation of linear elasticity is in material coordinates 
{X, f} , because the elastic properties of a continuous medium depend on 
the physico-chemical properties at every material point. I define linear 
elasticity as the theory one obtains when assuming (this is my version of 
Hooke's law) that, in material coordinates, 



s J j(X, f;*o) = c I 7 K I (X / fo)E K L (X / f;f ) , (55) 



20 I don't know of any text there this relation between an instanteous rotation velocity 
(v(t) and the associated finite rotation R(f; to) is given, other than my own Elements for 
Physics (Tarantola, 2006). 
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where the positive definite stiffness tensor has the symmetry 

ClJKL = cklij ■ (56) 

Note that, at any instant t , I write Hooke's law using the components of 
the stiffness tensor "frozen" at £o • 

In the laboratory coordinates, one has s'j = \ l j\hs l ] = \ l i?J jC 1 j k lEk l = 
A i I \JjC I j K L A k K \ L £ E k £ ,i.e., 

jj(x,t;to) = c i j k e {x,t)E k e {x,t;t ) , (57) 

where c) k e (t) = A'i(t) A' ; (£) \ k K (t) \ L e (t) c'jhik) . 

Appendix 11.3 analyzes the example of isotropic elasticity. It is there 
explained the well-known fact (e.g., Nowacki, 1986) that, while in the sym- 
metric theory, the isotropic stiffness tensor has two invariants, in the gen- 
eral theory it has three. 

In material coordinates, the Hooke's law implies the relation 

s'/pU^o) = c'/Vx^E/CX^o) , (58) 

but there is no simple relation between (t; to) and dj (£) . 

Needless to say, the theory here presented is just the mathematically 
simplest theory. Physical reality may suggest that the "constants" c 7 j K i,(fo) 
may, in fact, be functions of the temperature, the state of deformation (or 
the stress), etc. Also, one may need to replace the linear relation (55) by a 
more general relation. Using, for instance, some terms of a series expansion 
would lead to 

s J / (X,f;i ) = c I / K L (X,io)E K i (X,i;i ) 

+ d 1 j K l M n(X, to) E^ L {X,t;to) Em N (X, £;£q) + • • • , 

d , . . . being appropriately defined tensors. 



6 Geodesic movements 

We shall say that a movement is geodesic 11 if the material components of the 
deformation velocity tensor are not time-dependent, d^(X, t) = dj^(X) . 

21 In the symmetric theory, it also has the symmetries cjjki = C]ikl = c ijlk ■ 
22 We see in section 8 that such a movement actually corresponds to a geodesic path in the 
configuration space. 
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Then the matricant series (51) just becomes the exponential function, and 
one has 

A(X,t;t ) = expQ(X / f;f ) , (60) 

where Q(X,t;to) is the tensor whose components in the material coordi- 
nates are 

Q/(X,t;t ) = (t-t ) d/(X) . (61) 
The strain being E = log A , one then has E^(X,t;to) = Qi'(X, t; to) , i.e., 



E/(X,f;f ) = (t-t ) dJ{X) . (62) 



So, in a geodesic movement, the components of the strain are just proportional to 
the components of the deformation velocity. In particular, one has 

£j{X t t;h) = d/(X) . (63) 

Note that this identity between strain rate and deformation velocity is only 
valid for geodesic movements. 

In the laboratory coordinates, = A^A-'jE/ = (t — to) A 7 , A J j d^ , 

i.e., 

Ei'{x,t;t ) = (f-f ) dj{x,t) . (64) 

where dj{t) = A f f (f) A' 7 (f) . Note that even if the movement is geo- 
desic, the components of the deformation velocity in the laboratory coordi- 
nates are time-dependent. 

7 Power and energy 

The volumetric power produced by the causes of the motion is 

w{X,t) = s'jiXj) dJ{X,t) , (65) 
and the volumetric energy cumulated between instant to and instant t is 

u(X,f;f ) = [ dt'w(X,t') . (66) 
J t 

Let us first examine the case of geodesic movements (section 6). Then, 
using Hooke's law (55) and the geodesic strain relation (62), one arrives at 



w 



(X,t) = (t - to) c I j K L (X,t )d K L (X)d I J(X) . (67) 
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Therefore, 

«(*;*□) = { -^^c I j K L (to)d I J d K L , (68) 

a relation that, using again expression (62) for the geodesic strain, can be 
written 



u(t;t ) = icV^^E/^o) E/(^o) 



(69) 



So, for geodesic movements, the volumetric elastic energy is a quadratic function 
of the strain. 

The obvious question is: when the movement is not geodesic, does this 
elastic energy depend on the path of the movement in the configuration 
space? Should the answer be negative, then expression (69) would have 
general validity. This is an open question, whose answer will require to 
further extend the already known properties of the matricant. 



8 Configuration space 

The configuration space of a deforming elastic medium, identified with 
GL + (3), was introduced with detail in Tarantola (2006), where some pic- 
torial representations, similar to those in figures 2-4 are presented. 

Every (asymmetric) deformation tensor A (as introduced in section 5.4) 
corresponds to a point of GL + (3) , so a general movement A(f ) is an arbi- 
trary path in GL + (3) . 

Let me now present some basic notions on the geometry of the Lie 
group manifold GL(n) . It is well-known that Lie group manifolds have 
a connection and a metric. The connection is not symmetric, so it is equiva- 
lent to say that Lie group manifolds have a metric and a (Cartan's) torsion. 
In reality, the torsion of a Lie group manifold is totally antisymmetric, this 
implying that geodesic lines and autoparallel lines coincide. Therefore, one 
can limit oneself to talk about geodesies. It is also well-known that not all 
the points of GL(n) can be reached geodesically from the origin. 

A matrix M of GL(n), has "entries", say M a b . When choosing these 
n 2 entries as a coordinate system on GL(n) , there is the associated natural 
basis at every point. Tarantola(2006) demonstrates the following property: 
if M belongs to the part of GL(n) that is geodesically connected to the origin, 
then m = log M is real, and the entries of m are the components on the natural 
basis at the origin I of the group of the vector of the tangent space ("algebra") that 
is tangent to the geodesic connecting I to M (on the group) and whose norm is 
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Figure 2: In two-dimensional elasticity, the configuration space is a part of 
GL + (2), and, if the deformations are volume-preserving, a part of SL(2) . 
A (partial) section of SL(2) is represented here, together with some its 
geodesies (thin black lines). See Tarantola (2006) for details. Each point 
of SL(2) can be seen as a transformation: that transforming one vector ba- 
sis into another. Here, the basis at the bottom-left (origin of the group) is 
transformed into the other bases represented. The points in the yellow area 
can no be reached geodesically from the origin. The logarithm of the ma- 
trices in the yellow area are not real (the components of the logarithm are 
complex numbers). 
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Figure 3: In two-dimensional elasticity, the part of SL(2) that is geodesi- 
cally connected to the origin is made of the matrices whose logarithm is 
real. Let A be such a matrix. It can be identified to the deformation tensor 
of the text. Then, the strain E = log A is defined (it is real). The symmetric 
part of the strain represents a deformation, and the antisymmetric part, a 
(Cosserat) rotation. Both, the deformation and the rotation are represented 
in this figure. In the grey area, points that belong to SL(2) but not to the 
configuration space. Of course the configuration space contains all possible 
deformations and rotations. 
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Figure 4: Continued from figure 3. The curved line represents an arbitrary 
movement. The movement velocity A (denoted "deformation velocity ten- 
sor" in the figure), is the tangent to any point of this line (in reality after its 
transport to the origin of the group). In the terminology of Tarantola (2006), 
this is the declinative of the path. The final configuration of the movement 
is the big pink dot. The strain associated to this final configuration is the 
dashed geodesic segment (this is the logarithm of the A representing the 
point). For the points of the configuration space not geodesically connected 
to the origin the strain is not defined (and the logarithm would not be real): 
physically, these configurations can not be accessed from the reference con- 
figuration by a stress of the form A s , where s is a fixed ("nominal stress") 
and A some parameter (for instance, the time t ) varying from zero to in- 
finity. 
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equal to the length of the geodesic. From where the logarithmic definition of 
the strain. 

Now, consider a point A of GL + (3), that belongs to the configuration 
space 23 . Then, E = log A is real, and can be interpreted as the (dashed) 
geodesic segment in figure 4. The decomposition of the strain into its sym- 
metric and its antisymmetric part, 

E = E + E , (70) 

allows to consider the (symmetric) deformation tensor 

D = expE (71) 
— that is a "standard" deformation — and the (orthogonal) rotation tensor 

R = expE (72) 

that corresponds to the Cosserat micro-rotations. This decomposition of A 
into a (symmetric) deformation D and an (orthogonal) rotation R is at the 
basis of the representation of the configuration space in figures 3 and 4. 
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11 Appendices 

11.1 The matricant 

What follows is an exposition of the notion of matricant, as exposed by 
Gantmacher (2000). I limit myself to small adaptations of notations and 
of language. 

Letting X(f) and P(t) be two time-dependent matrices, Gantmacher 
considers the differential matrix equation 

X(t)Xit)- 1 = P(t) . (73) 

Here, P(t) is assumed to be "a continuous matrix function of the argument 
t in some interval (a, b) ". A solution to the system (73) is sought such that 
for some to in the interval (a, b) , the solution satisfies X(to) = I. Such 
a solution is determined by "the method of successive approximations". 
The successive approximations X^{t) (k = 0,1,2, ... ) are found from the 
recurrence relations 

M^Xt-iit)- 1 = P(f) (k = l,2...) , (74) 

where Xo(t) is taken equal to the identity matrix I. Setting K^(fo) = 
I (k = 0, 1, 2, . . . ) , one may represent X^{t) in the form 

X k {t) = 1+ f dt' P(f') X^O , (75) 

thus, X (f) = J, Xi(f) = 1+ f t [dt' P(t'), X 2 (t) = 1 + Jl Q df P(f)(l + 

fl dt" P(t")) = I + J t [ dt' P(f ) + H dt' P{t>) f< dt" P(t") , etc. Then, 
Gantmacher proves that this series, 

t t t' 

n(f;f ) = 1+ f dt'P(t')+ f dt' P(f') f dt"P(t") + ... , (76) 

J to Jtg J t(, 

is absolutely and uniformly convergent in every closed subinterval of the 
interval (a, b) , and, therefore, constitutes the solution of (73). 
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That the sum (76) is the solution of (73) is verified by a term-by-term 
differentiation. "This term-by-term differentiation is permissible, because 
the series obtained after differentiation differs from (76) by the factor P(t) 
and, therefore, like (76), is uniformly convergent in every closed interval 
contained in (a, b) " . As already anticipated in expression (76), this "nor- 
mal" solution (often called the matricant) is denoted Q(f;fo) ■ Gantmacher 
explains that every other solution is of the form 

X(t) = O(f;f )C , (77) 

where C is an arbitrary constant matrix. Gantmacher says that it follows 
from this formula that every solution, in particular, the normalized one, is 
uniquely determined by its value for t = to. 

The representation of the matricant in the form of such a series was first 
obtained by Giuseppe Peano (Peano, 1888). The matricant theory is used 
in seismology, typically for the propagation of wave fields in depth: I first 
learned about the matricant when reading Brian Kennett's book (Kennett, 
1983). 

Property #1: One has Q(f;f ) = C\(t;h) Q(h;t ) . 

Property #2: One has Q(t;t )(P + Q) = Q(t;t )(P) Q(t;t )(S) , with 
S = Cl{t ;t){P) QCl{t;t ){P). 

Property #3: One has det Q(f; f ) = exp ( dt' tr P(t')) . 
Property #4: If P is constant, Q(f;f ) = exp ((f - f ) P) . 

11.2 Transpose and adjoint 

Let V be a linear space, *V its dual. The mathematical definition of the 
dual of a linear space is abstract 24 , but we only need here the most basic 
of its properties: if V is a space of vectors with components v 1 (in some 
vector basis), then X = * V is a linear space of objects (forms) with com- 
ponents Xi (in some form basis), so that the expression 

(X, V) EE X,V EE £>V (78) 

i 

makes sense. This is called the duality product. 
A relation like 

w = Sv ; w a = S a i v i (79) 
24 It is the linear space containing all the linear forms over V . 
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defines a linear mapping from a linear space V into a linear space W . The 
transpose of the mapping, denoted f S , is, by definition the linear mapping 
from ¥ = * W , the dual of W , into X = * V , the dual of V , such that 
the relation 

( f Sy,v) = (y,Sv) (80) 
holds in general. Explicitly, this is 

( f Sy)^ = y«(Sv)« . (81) 

While the components of S were denoted using the indices S a , it is con- 
venient to denote f S using the same symbol S , and just changing the po- 
sitions of the indices: Sj a . The condition in equation (81) then becomes 

(Si a y a ) v { = y a (S" t v { ) . (82) 

It is clear that for this relation have gereral validity, one must have 

Si a = S a i , (83) 

that is the relation holding between the components of a linear mapping 
and its transpose. Practically, excepted for a "replacement of the indices" there 
is no difference between a mapping and its transpose. 

In the same context, assume now that the two spaces V and W are, 
in fact, scalar product vector spaces, i.e., assume that that there exist two 
metric tensors gjj and defining the two scalar products 

( vi , v 2 ) = gij v\ v[ ; ( wi , w 2 ) = 7 a/J w\ w\ . (84) 

Then, in addition to the transpose, one can introduce the adjoint, denoted 
* S , that is, by definition the linear mapping from W , into V , such that the 
relation 

(*Sw, v) = ( w, Sv) (85) 
holds in general. Explicitly, this is 

# ; -(*Sw)V = lu ^{Svf = 7«^ a (S^^) , (86) 

i.e., using the relation (82) involving the transpose, 

£ !7 fSw)<V = (S/t^O^ . (87) 

In order for this to hold unconditionally, one must have 

(*Sw) ! ' = *S a 'w a , (88) 
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with 

*S a > = -y^Sfg' 1 . (89) 

The situation found in the text is a special case of this, where the map- 
ping S is an endomorphism (mapping a linear space into itself), so there is 
only one metric. 

11.3 Fourth-rank isotropic (asymmetric) tensors 

Here, the notion of fourth-rank isotropic tensor is discussed, without any 
particular reference to elasticity or viscosity. 

The viscuous of elastic invariants (eigenvalues of the fourth-rank isotro- 
pic tensor) will typically be a function of the material point. If using the 
material coordinates, we shall then face the functions 

A K = T K (X) ; A F = r„(X) ; k w = r w (X) , (90) 

while, if using the laboratory coordinates, we shall face the functions 

K = y K (x,t) ; A,, = j }l (x,t) ; A w = y u (x,t) , (91) 

related to the previous ones via 

A K = J K (x,t) = r K (o(x,f)) 

\ v = 7p (x,t) = r H (0(x,f)) (92) 
Ao; = 7w(x,t) = T w (®(x,t)) . 

Note the time-dependency of the functions in the laboratory coordinates. 

In material coordinates, the components of a fourth rank isotropic ten- 
sor are 

c l j\{x,t) = i K {x)k l j\ 

+T H (X)m I J K L (X,t) (93) 
+T w (X)a I J K L (X,t) , 

with the three orthogonal projectors 

UK _ 1 :I :K 

J L ~ 3<> / L 

m l j\(X,t) = l(g IK (X,t)gj L (X,t) + S I L Sj K )-lS I J S K L (94) 

a l f L {X,t) = \{g lK (X,t)g ]L {X,t)-5 l L 6f) . 



30 



In laboratory coordinates, the components of a fourth rank isotropic 
tensor are 

<?j k t{x,t) = y K (x,t)iej k i 

+ 1}l (x,t)m i j k i (x) (95) 
+7w(x,t)a'i k e (x) , 

with the three orthogonal projectors 

k l j k e = \ S l j S k £ 

= \{g ik {x)g j ,{x)+8 i ,6 j k )-\8 i j 5 k l (96) 
«?<(*) = l(g lk (x) gj( ix)-5' e 5 j k ) . 

Note that while in the material coordinates, the time-dependencies are 
in the components of the metric, gij(X,t) and g^(X,t), in the material 
coordinates they are in the functions j K (x,t) , jp (x, t ) , and ju>(x,t) . 

These components are related via 

c l ] K L = \ ! i tij \ K k X l h eft ; df t = A'j \ J j A k K A L , c l ] K L , (97) 
as it should. 



See the next two pages for the Tables of Formulas. 
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11.4 Tables of formulas 
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material coordinates {X,t} 
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motion: 
coefficients: 
metric: 
material velocity: 
deform, velocity: 
property #1 : 
[vorticity] : 
viscosity: 


deform, gradient: 
sq. def. tensor: 
property #2: 
deformation tensor: 
strain: 
elasticity: 


property #3: 
property #4: 
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laboratory coordinates {x, t} 
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material coordinates {X, f} 
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motion: 
coefficients: 
metric: 
material velocity: 
deform, velocity: 
property #1: 
[vorticity]: 
rotation velocity: 
asym. def. velocity: 
viscosity: 


asym. deformation: 

asym. strain: 
elasticity: 


property #2: 
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